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Abstract 

q ' We study relativistic mean-field models with hadron masses and coupling constants 

depending self-consistently on a scalar meson field. We demonstrate that by the field 
redefinition some models can be equivalently transformed into each other. Thereby 
the large variety of scaling functions for masses and couplings can be reduced to 
^ ■ a restricted set of functions with a constrained dependence on a scalar field. We 

show how by choosing properly the latter scaling functions one may stiffen or soften 
the equation of state at high densities and simultaneously increase the threshold 
density for the direct Urea process without any change of the description of nuclear 
matter close to the saturation density. The stiffening of the equation of state might 
be motivated by recent neutron star mass measurements, whereas the increase of 
the threshold density for the direct Urea process is required by the analysis of 
neutron star cooling data. We demonstrate that if a rho meson is included in a 
mean-field model as a non-Abelian gauge boson, then there is a possibility for a 
charged rho-meson condensation in dense nuclear matter. We show that such a 
novel phase can be realized in neutron star interiors already for sufficiently low 
densities, typically ~ 3-^4 no, where uq is the nuclear saturation density. In the 
framework of the relativistic mean field model the new phase arises in a second- 
order phase transition. The appearance of a p~ condensate significantly alters the 
proton fraction in a neutron star but changes moderately the equation of state. The 
neutrino emissivity of the processes involving a p~ meson condensate is estimated. 
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1 Introduction 



The baryon matter at densities relevant for neutron star (NS) interiors is con- 
venient to describe in terms of relativistic mean-field (RMF) models, cf. [1] 
and references therein. In these models the nucleon mass decreases with an 
increase of the nucleon density due to coupling to a scalar field, whereas the 
effective masses of scalar (a), vector (u>) and isovector (p) mesons, and their 
coupling constants, are usually assumed to be density and field independent 
quantities. The RMF models provide a simple and economical tool for con- 
struction of the equation of state (EoS) for dense baryonic matter. By tuning 
parameters of a RMF model at a saturation density one may get an EoS sim- 
ilar to the one that follows from more involved microscopic calculations, e.g., 
as that of Urbana-Argonne group (Al8+5f+UIX*) [2]. The latter calculation 
provides probably the most realistic EoS for nuclear densities up to 4n , but 
it uses a non-relativistic potential and violates causality at higher densities. 
As an extrapolation to higher densities RMF models are practical. However 
all RMF models have the very same an unpleasant feature. They produce a 
large fraction of protons in the NS matter, cf. [3,4]. This permits a very effi- 
cient cooling process, direct Urea (DU) process n — > peu, for n > n™, with 
a low threshold density 

^crit ^ ^ ^0* 

Ref. [5] emphasized this fact as an ob- 
vious shortcoming of RMF models. Oppositely, ref. [6] and some subsequent 
works [7] tried to use this fact to develop a so-called DU-based scenario of a 
NS cooling (based on ~ (2 4- 3)n ). 

The NS cooling " data" (surface temperatures of NSs as a function of the NS 
age) can be subdivided in three groups, corresponding to "slow cooling", "in- 
termediate cooling" and "rapid cooling" of neutron stars. If the density in 
the NS central region exceeds n™, which corresponds to the mass M^Y, the 
star cools down very fast, radiating mostly from its central region. Thus, the 
NSs with M < M°;Y cool slowly, whereas the stars with the mass M only 
slightly above would be already characterized by a very rapid cooling. 

Such a scenario was criticized in [8,9,10]. It is doubtful that many NSs, be- 
longing to an intermediate cooling group, have the very similar masses. It 
has been already experimentally established that NSs may have essentially 
different masses. E.g., masses of two NSs in double NS systems are mea- 
sured with the very high precision, MB1913+16 = 1.4408 ± O.OOO3M , and 
Mjo737-3039B = 1 .250 ± O.OO5M [12]. 

Authors of many RMF models do not care about the DU threshold at all 
and obtain so low values of that < 1 M & . However according to 

the standard mechanisms only the NSs with M > 1 M & can be formed in su- 
pernova explosions. The mentioned above microscopic Urbana-Argonne model 
produces nj? r ^ ~ 5 no, that corresponds to a heavy NS with ~ 2M Q . 

If one used a standard RMF model fitted in [11] to reproduce the Urbana- 
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Argonne EoS for n < 4uq, but, as well as other RMF models, producing a 
low DU threshold (n££ ~ (2.6 -=- 2.7) no), one would conclude that the ma- 
jority of the NSs seen in soft X rays are low mass objects with the mass 
M < 1.3 Mq. The latter value is however below the averaged value of the NS 
mass (M ~ 1.35 ± 0.04 M ) measured in NS binaries [13]. 

Another problem is that there might exist heavy NSs. Recent measurements 
of NS masses in binary compact systems yielded Mj 75i+iso7 = 2.2 ± O.2M 
[12] and M 4U i700-37 = 2.44 ± 0.27 M [14]. The latter object might be a 
black hole due to a lack of pulsation but the former object, J0751 + 1807, is 
definitely a NS. The confidence level of 95% still allows to have a smaller 
limiting NS mass than 2 M (however the limiting NS mass is surely higher 
than 1.6 Mq). Additional information on heavy NSs may be obtained from 
quasi-periodic oscillation (QPO) sources. These are NSs that are X-ray sources 
at frequencies of the orbiting accreting matter. If QPOs originate from the 
innermost stable orbit [15] of the accreting matter, their observed values imply 
that the accreting NS has the mass ~ 2.2 M in the case of 4U1820-30. The 
value 2.2 M is the limiting mass of a NS following the Urbana-Argonne EoS. 
Unfortunately this EoS is acausal at high densities. Ref. [16] suggested an 
analytic fit of the Urbana-Argonne EoS for n < 4n , which respects causality 
for larger densities. We will denote this parameterization as the HHJ EoS. 
As a price, the limiting mass decreased down to ~ 2M . Moreover, all phase 
transitions, like pion and kaon condensations, and quark matter, which are 
assumed to be possible in NS interiors, may essentially soften EoS decreasing 
further the limiting mass of a NS 1 . Thus modern microscopically motivated 
EoS and possible phase transitions in NS interiors depend crucially on one 
precision measurement of a heavy NS with M > (2 4- 2.2) M . 

Within a standard RMF model it would be hard to describe heavy NSs with 
M > 2.2 Mq without either an essential increase of the compressibility param- 
eter K or a decrease of the effective nucleon mass at the saturation density. 
Both values are rather constrained by the atomic nucleus data. The problem 
originates from the fact that the RMF model parameters are fixed at the sat- 
uration density and no free parameters remain to control the pressure and the 
symmetry energy at higher densities. The problem could be in principle solved 
by introducing new terms with new fitting parameters. 

There are strong indications that meson masses as well as the nucleon mass 
depend on the nucleon density. A general trend for such a dependence is 
governed by the property that masses of all hadrons except Goldstone bosons, 



1 The Urbana-Argonne calculation [2] already includes neutral pion condensate for 
n > 1.3 no that only slightly softens EoS in their case. 
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like pions and kaons, 2 should decrease, when the nucleon density n increases 
towards a critical point of the chiral symmetry restoration [18]. According to 
the conjecture of Brown and Rho [19] the nucleon mass and the masses of 
vector uj , p and scalar a mesons may obey an approximately universal scaling 
law (a so-called "Brown-Rho scaling") 

m* m* m* m* N 
m p m^j m CT m N 

where asterisks indicate in-medium masses. According to [19] <&br decreases 
with the increase of the density. In ref. [20] an attempt was made to in- 
clude density dependent masses in the RMF model in order to incorporate 
the Brown-Rho scaling, explicitly. A theoretical disadvantage of this model is 
that the density dependence, $br — (1 + yn/n )~ 1 , y ~ 0.28, was inserted by 
hands into the model Lagrangian rather than it followed from the correspond- 
ing equations of motion. Insertion of a density dependence into parameters of 
a mean-field Lagrangian violates conservation of the energy-momentum and 
spoils the thermodynamic consistency. Therefore, such models require a spe- 
cial treatment [21]. Ref. [20] also realized that the scaling of meson masses (1) 
would result in a dramatic stiffening of the EoS, if one does not simultaneously 
scale the corresponding coupling constants in a similar way. As an argument 
for such a procedure, authors relied on the description of the flow in heavy 
ion collisions. Below we will show that this prescription relies on equivalence 
of some RMF models that allows to use different interpolating mean fields 
without changing physical quantities on a mean-field level in homogeneous 
matter. 

In ref. [22] the density dependence of effective hadronic masses was incor- 
porated into the mean-field Lagrangian via extra terms describing a hadron 
interaction with a scalar field. It resulted, however, in an increase of the p and 
uj meson masses with increase of the nucleon density in disagreement with (1). 
At the same time the a mass was assumed to be constant. Thus, the inclusion 
of the hadron mass scaling into RMF models deserves a further considera- 
tion. Some other works, cf. [23], played with a field-dependent parameters of 
generalized mean field models, permitting changes of both masses and cou- 
plings. Equivalence between different parameter choices was not realized in 
these papers. 

Appropriate knowledge of the EoS of hot, dense, and isospin asymmetrical 
nuclear matter is essential for understanding the results of heavy-ion collisions. 
Flow, particle yields, fragmentation may give some constraints on the EoS and 
the symmetry energy, cf. [24,25]. Often these constraints disagree with those 

2 The kaon effective mass may also decrease with the density, whereas the nature 
of this decrease is different [17,1,11]. 
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obtained from NS physics, e.g., yielding very low DU threshold density. The 
proper choice of the EoS is particularly important for the CBM (compressed 
baryon matter) experimental program at the future facility at GSI. Precision 
measurements of heavy ion collisions are planned in the whole energy range 
from ~ 1 GeV/A up to ~ 40 GeV/A. 

To describe correctly the symmetry energy of the nucleon matter one needs 
to include the p-meson field into a mean-field Lagrangian. The p-meson field 
is traditionally assumed to have only one mean-field component, (p M ) a = 

<Vo<5 a3 Po 3 \ where \x = 0,1,2,3 is the Lorentz index, and a = 1,2,3 is the 
isospin index. This field produces then the expected contribution to the energy 
density oc (n p — n n ) 2 , where n p and Ti n £1X6 the proton and neutron densities. 
In ref. [26] it was realized that, if the p-meson is introduced as a non-Abelian 
gauge boson, e.g., in the framework of the hidden local symmetry model [27], 
then at a density n > n p c the charged p-meson condensation occurs. In a new 
phase the symmetry energy grows as oc \n p — n n \ only. The critical density 
was found to be n£ ~ ^p 3 , being sensitively dependent on the value of the 
in-medium mass of the p-meson (m*). The phenomenon of the charged p- 
meson condensation is similar to the condensation of gauge bosons in QCD [28] 
(gluon condensate) and to electro-weak [29] (I¥-boson condensate) sectors of 
the Standard model. 

In this paper we first construct a generalized RMF model with all parameters, 
being field dependent (sees. 2,3). We derive the energy density functional, and 
show how various standard cases can be reproduced. We discuss the equiva- 
lence relations among some cases. Then we demonstrate shortcomings of the 
models with a hadron mass scaling without a scaling of coupling constants. 
We show that the RMF models with the universal scaling for masses and 
coupling constants can be reduced to the models without any scaling. Such 
models do not solve, however, the problem of the low n^. Then we allow 
for a different scaling of hadron masses and coupling constants in order to 
generate a model that produces a stiffer EoS at high densities. This model 
supports a large limiting mass of the NS and rather high critical density for 
the DU process. On the other hand, it uses appropriate values of compress- 
ibility and effective nucleon mass at the saturation nuclear density. Then we 
demonstrate another example of the RMF model with a non-universal scaling 
of masses and coupling constants that fits well the HHJ EoS including the 
correct density dependence of the proton concentration. In sec. 4 we investi- 
gate a possibility for a charged p-meson condensation. We demonstrate that in 
the framework of our generalized RMF models the charged p-meson conden- 
sation occurs as a second-order phase transition. Emissivity of the NS cooling 
processes involving the charged rho condensate is estimated in sec. 5. Results 
are summarized in the concluding section 6. In Appendix we show how one 
can introduce inhomogeneous electric potential. Throughout the paper we use 
units h — c—1. 
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2 Generalized relativistic mean-field model 



2. 1 Mean field Lagrangian 



We start with a generalized RMF model, where effective hadron masses and 
coupling constants are assumed to be field-dependent from the very beginning. 
The Lagrangian density is given by 



C — Cn + Cm + C\ , (2) 
C N = a N ^ N (iD ■ 7) - rn N <j) N * N ^ N , 

D tl = d ll + ig u ,XuU ll + -j9 P XpPixT 1 (3) 
C M = a a <p a — U{a) 

= - 9^ , = d u p p - d^p v + g' p x' p [p P x p,,] , (5) 
C l = J2^i[Kl-d)-m l ]^ l . (6) 

Here = (^n,^ P ) T is the isospin doublet of the nucleon bispinors, a, u; M , 
p M are the a, a; and p meson fields, r is the isospin Pauli matrix. The part Ci 
describes the contribution of the light leptons (electrons and muons, I — e,p). 
For the sake of simplicity we will neglect hyperons in the Lagrangian (2), which 
can contribute, when the hyperon Fermi seas are filled in the NS matter. 
An extension for hyperons is straightforward, but is associated with large 
uncertainties due to poorly-known coupling constants of hyperons to mean- 
fields, cf. [1,11]. The meson sector in (2) is treated on a mean-field level as 
described in Ref. [1]. 

The scaling parameters of the model, the "dielectric" constants a^, renormal- 
izations of the coupling constants and masses, %i an d <f>i, are dimensionless 
functions of the fields and effective couplings. The potential U(cr) represents 
the possible self- interaction of the a field, which was suggested in ref. [3]. 
In the general case the scaling parameters can depend also on (uj fl uj fl ) and 
(PfiP^ 1 ) yielding corresponding terms in the potential. For instance, writing 
4>t = 'fitW) + C^m^V^ w e recover the (uo p uo^) 2 terms proposed in ref. [30] to 
improve the description of finite nuclei within the RMF model approach. A sys- 
tematic construction of mean-field Lagrangians with high-order terms based 
on constraints from the chiral symmetry of QCD can be found in ref. [31]. In 
our treatment we will neglect a possible dependence of the potential and the 
parameters on the uj and p meson fields, elaborating only the a-field depen- 
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dence. We assume that the density dependent scaling functions a at and 4>n 
are basically due to a nucleon-sigma-meson interaction. Thus, we can write 



a N = a N (g a Xa(?) , 4>n = <j>N(.9<,X<r°) ■ ( 7 ) 

Simplifying we will also assume that all other scaling functions, x^,pi &u,p,oi 
and the potential U are due to the same nucleon-sigma-meson interaction, i.e. 
are the functions of g a x a o. For the rho meson included according to the hidden 
local symmetry principles [27] one has g p = g . Despite this the scaling of g p 
and g may be quite different in dense nuclear matter, xp 7^ Xp- The scaling 
factor x.p arises due to the renormalization of the nucleon-sigma interaction in 
medium, whereas the factor x p is due to a non-Abelian interaction between 
fields and does not depend directly on the nucleon-field source. Thus within 
above assumptions we may take x p ~ 1. 

The nucleon density is given by the zero component of the Noether current 
n N = a N < tf N y N > . (8) 

2.2 Field redefinition in RMF models 



Describing an extended homogeneous system like a NS we may assume that 
mean fields of mesons do not depend on coordinates. Then the physical results 
obtained on a mean-field level would not change if all particle fields are rescaled 
by constant factors. For instance we may change meson and nucleon fields as 



-> ^n/ y/aN , cr -> a I y/bZ , Up -> uj ^faZ , p p -> pj ^/a~ p (9) 

and obtain the new nucleon and meson Lagrangians equivalent to those in 
(3,4) in a mean-field approximation 



C N = V N (iD- 7) - m* N ^ N 

Dp = d fl + ig u ,XuU fl + -g P X P PpT, (10) 



■>p,v w a "u, w'/jW _ P pu P m p P pP / , , \ 

4 + 2 4 + 2 ' 1 } 

ujpu = d u u^ - d p u v , p pu = d v p^ - d p p u + g p x p [pp x p v \ . (12) 
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We denoted U{o) = U(xa&)- The effective masses of particles introduced in 
this Lagrangian are equal to 



m* N /m N = (p N (xav)/a N (x*v) = ®n(x*°) , 




= <T,U,p, 



(13) 



where the dimensionless functions $at and depend on the scalar field in the 
combination x<t°- The scaling functions of coupling constants g a , g^, g p , g p are 
given by 



All quantities are now expressed in terms of the new renormalized a field. 

The thermodynamic potential density Q = —P, where P is the pressure, is 
related to the energy density by the standard equation 

BE 

E[n n ,n p ,ni;f,p l ,} = Y^tJ>ini + tt, P = ^P' i n i -E 1 ^ = — . (15) 

i i 1 

Chemical potentials enter the Green functions in the standard gauge combi- 
nations £j + p,i. Chemical potentials of the charged boson mean fields are also 
introduced by the gauge replacements — > Si + pi. 

For charged rho-meson fields above replacement implies 

P^v -> Ptiu + bp,,/, , Apfiu = p P ch 5uo[n 3 x p M ] - p p ch S p0 [n 3 x p v \ , (16) 

where {n 3 ) a = 5 a3 is the unit vector in the isospin space, a = 1, 2, 3 , and p p ch 
is the chemical potential of rho mesons. The explicit expression for Q is given 
in Appendix. 

After calculation of the nucleon contribution the energy-density functional 
takes the form 

E[n n , n p , m\ f, uj , fa] = E N [n n , n p ; f] + Ei[n e , n p ] 




(14) 



+EJ 

Tlni Tip ; f,u ] + E p [ 

^ni Hp, /, Pp] , 



(17) 



E N [n n ,n p ] f] 




Va(f) + U(f) 



2CI 




(18) 
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PF 



/(\ppP j 
-^-ymf + p 2 , 



i=e,fi q 



C 2 (n n + n p ) 2 
2m 2 N rjM) 

2c 2 



9u Xuj - 

(3)n 



Clin p + n n ) 



E p [n n , n p ; f, pp] — E Pjneut [n n , n p ; f, p ] + E p ^\n n , n p ; f, pp] , 



Ep >ncu t[n n , n p , f, Pq 



(3)i 



C 2 (n n - n p y 



8m 2 N r] p (f) 
m 2 NV P (f) 



2C 2 



9 P XpPo 



1m 2 N ri P (f) 



EpAnm n p , f, pp] = (p^ p» ~ P ( o ] P ( o ] 



We have introduced here the dimensionless field 



/ = 9aXaCr/m N , 



(19) 



(20) 



(21) 



(22) 



and the dimensionless coupling constants Cj = gi jmi , for % = a, uj p . The 
nucleon and lepton Fermi momenta pp,i are related to the corresponding partial 
density as n« = p%,J (3 for % = n,p,e, p . The contributions of mesonic fields 
entered with the scaling factors 



!*(/) =*?(/)/*?(/) , i = a,u,p. 



(23) 



Then one may use equations of motion for u and p- meson fields p M : 



d d 

tt[n n , n p ; f, uj ] = ^ — E[n n , n p ; f, u ] = , 



OUJq OOJq 

d d 
— Q[n n , n p ; f, pp] = — E[n n , n p ; f, p M ] = . 

OPp Opp 



(24) 



The first line in (24) yields 



Eu[n n ,n p ] f] 



C 2 (n n + n p ) 
2m 2 N r] LU (f) 



(25) 
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Traditionally one considers only a neutral p-meson field p^ disregarding a 
possibility for the presence of charged p-meson mean fields. Then the second 
line in (24) yields 



pf =g P X P {n p -n n )/{2my p ). (26) 
The corresponding contribution to the energy density is 

77 r n C 2 {n n - n p ) 2 

E P [n n ,n p J]= 8mjU(/) . (27) 

As one can see, the source of the isovector p-meson field Pq is the isovector 
nucleon density (n p — n n ) . Therefore for the isospin symmetrical matter (n p = 
n n ) we always have E p — . 

The main merit of this section is in the following. Eqs. (17-27) demonstrate 
that instead of twelve scaling functions of a scalar field, a,N,cr,u>,p, Xp,w,a, <pN,au>,p 
and U(cr), which enter the initial Lagrangian (2), the energy density func- 
tional depends actually only on four particular combinations of these func- 
tions, r] a ^ iUJ (f) and U(f). Note that the dependence on the scaling function 
r\ c can be always presented as a part of the potential U and absorbed in it 

by the replacement U — > U + 2 N c2 (1 — Va(f)) > an d vise versa the potential 
U can be absorbed in r] a . Thus only three independent functions enter the 
energy-density functional. The expression (17) demonstrates explicitly equiv- 
alence of mean-field Lagrangians for constant fields with various a, x an d <fi 
parameters if they correspond to the same functions r) PjU (f) and r\ a (either 
U(f)). In sec. 3 we discuss possible constraints on the / dependence of the 
functions. We would like to stress that the field / relates now to a scalar field 
a in a non-trivial way determined by eq. (22). 

In order to obtain final expression for the energy of the system, one should 
still use equation of motion for the scalar field /, 



d d 

— n[n n ,n p \f\ = —E[n n ,n p \f\=Q. (28) 



2. 3 Determination of parameters of a RMF model 



Parameters of the RMF model, C a , C u , C p , and parameters of the potential U 
are to be adjusted to reproduce the nuclear matter properties at the saturation 
and have rather broad uncertainties, cf. [32,33]. One usually takes the satura- 
tion density n — 0.16 ±0.015 fm~ 3 , the binding energy per particle, ee(no) — 
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-15.6 ± 0.6 MeV, the compressibility modulus K(n ) ~ 240 ± 40 MeV, the 
effective nucleon mass m* N {no) /mjv — 0.75 ± 0.1, and the symmetry en- 
ergy coefficient a sym (n ) = 32 ± 4 MeV. Some models exploit parameters 
beyond even these large error bars. E.g., the original Walecka model [34] used 
n ~ 0.191 fur 3 producing K ~ 540 MeV and m* N /m N ~ 0.56. 

To facilitate the comparison among different RMF models which we study 
below we will use the same basic input parameters 

n = 0.16fnr 3 , e B = -16 MeV, a sym (n ) = 32 MeV, 

m N = 938 MeV . (29) 



The saturation density and the binding energy are related as 



dE[n; f] 



dn 



= — E[n ; /(n )] = m N + e B , 

"o,/Oo) n ° 



(30) 



and the compressibility modulus is given by 



K = 9n 



d 2 E 




/ d 2 E 


)' 


~d 2 E 




-1" 


dn 2 


"0, /(«(>) 


\dndf 


no,f(n )/ 


df 2 


n J (no) _ 





(31) 



Here f(n ) is a solution of eq. (28) at the density n p = n n = n /2. The 
coupling constant C p is determined from the symmetry energy coefficient of 
the nuclear matter 



, . n d 2 _, r . C 2 n 7r 2 n . , 

a sym (n) = ^-^E[n - n p , n p \ = — ^— + , (32) 

Pfn is the nucleon Fermi momentum in isospin-symmetric matter. 

Generic Lagrangians, and £^ in (2) and the energy density functional 
(17) reproduce various types of mean-field models. Consider several examples. 
The choice 



<lN,(T,u),p — XN,a,LU,p — 4><r,u),p — 1 , 4>N ~ 1 ~ Qa & /m N , (33) 

x P = u = o, 

yielding the energy-density functional (17,18,25,27) with r) c ^,p = 1, and $at = 
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1 — /, corresponds to the standard Walecka model (W) 3 with the minimal 
number of free parameters (since U — 0) [34]. For (29) we find the parameter 
set 

W: Cl = 329. 70, C 2 U = 249.40, C 2 p = 68.09 , (34) 
which produces 

K ~ 553 MeV, m* N (no)/m N ~ 0.54. (35) 

Thus this minimal model does not allow to fit appropriately the values of the 
nucleon effective mass and compressibility modulus. In order to cure these 
problems, it was suggested in ref. [3] to introduce a scalar-field self-interaction 

U(f) = m%(b f 3 /3 + c / 4 /4) = I (m N g a <x) 3 + | (g a a) 4 . (36) 

Two additional parameters, b and c, make the model (called the modified 
Walecka model (MW)) able to accommodate realistic values of the nuclear 
compressibility and effective nucleon mass. An extra attention should be paid 
to the fact that the coefficient "c" should be positive, otherwise there is no 
stable vacuum for the Lagrangian (4,11) 4 . The solid line in Fig. 1 shows 
values of K and m* N (n ), for the boundary case c = 0. Above this line c > 0, 
and below c < 0. With a positive c, low values of m* N (rio) correspond to 
unrealistically high values of the compressibility K. 

Using as input (29) and 

K = 270 MeV, m* N (n )/m N = 0.8 (37) 
we obtain 

MW : Cl = 189.94 , C 2 = 90.768 , C 2 p = 100.18 , 

b = 6.3714 x HT 3 , c = 1.6288 x 10~ 2 . (38) 

The choice, cf. (2), 

3 Inclusion of the neutral field generalizes the model to describe asymmetrical 
nuclear matter. We call the models, that include p-meson fields, by the same names 
as those without p-meson fields. 

4 Some parameter sets used, e.g., in [1], do not respect this constraint dealing thus 
with a metastable state. One then assumes that a mean field energy is valid only 
near a local minimum. 
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m > ) /m N 

Fig. 1. The line shows the compressibility K of the isospin symmetrical matter vs. 
the effective nucleon mass m* N (no) corresponding to c = for the MW model, input 
set (29). 

0>N = 1 + 9a^/ m N , a a,u,p = 0JV,cr,aj,p = XN,(T,u,p = 1) (39) 

x P = u = o 

corresponds to the original Zimanyi-Moszkowski model with the derivative 
coupling (ZM) [36]. In eqs. (17,18,25,27) we have then U = 0, r\ a ^,p = 1, 
^ct,lu,p = 1 and $at = 1/(1 + /) • For the ZM model with the input parameters 
(29) we obtain 

ZM : Cl = 169.36 , C£ = 59.055, C 2 p = 104.56 . (40) 

This model reproduces realistic values of the compressibility modulus and the 
effective nucleon mass at n , 

K ~ 224.5 MeV, m* N (n )/m N ~ 0.85 , (41) 

without introduction of extra free parameters. As the model, which eliminates 
successfully the deficiencies of the original Walecka model keeping the minimal 
number of parameters, the ZM model has gained a lot of attention and has 
been widely used. Notice that, the ZM model can be also considered as a 
generalization of the W model done with the introduction of a scaling of the 
<jN coupling 

9*a = 9aXa = 9a/ (1 + g^/m N ) . (42) 

To close this section we note that eqs. (13) and (14) demonstrate that for con- 
stant mean fields under consideration "dielectric constants" a,i can be always 
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absorbed into effective masses and effective coupling constants. Therefore, be- 
low we will put cij = 1 and deal with the functions $j and x% only. 



3 RMF models with dropping hadron masses 

In this section we construct RMF models which support a decrease of hadron 
masses with an increase of the baryon density Simplifying we follow the orig- 
inal Brown-Rho conjecture [19] about the hadron mass scaling, which in our 
case imposes 



Thus for the generalized W and MW models we would have $ = = 
$mw(/) = 1 - / and for the generalized ZM model, $ = $zm(/) = 1/(1 + /) ■ 
We recall that / is related to the a field via (22). 

We have to note that due to non-linear interactions effective masses obtained 
on a mean-field level should not be in general identified with positions of 
the poles in propagators of corresponding excitations. 5 Following Brown and 
Rho [19] we apply a mass scaling on a mean-field level. 

3. 1 RMF models with decreasing hadron masses and fixed coupling constants 

Assuming that coupling constants do not change in the medium, i.e. x% — 1) 
we obtain the following scaling factors in the energy- density functional 



For the generalized W and MW models we have $ = $w = $mw(/) — 1 — / 
and for the generalized ZM model, <E> = $zm(/) = 1/(1 + /)■ In terms of 
a a field these functions have familiar forms $mw(°") — 1 — 9a°'/ m N and 
$zm(c) = 1/(1 + gaO-/m N ) . 

The resulting functional possesses several unpleasant features. First, effective 
hadron masses do not decrease monotonously with the density increase. De- 
creasing at small densities, they start to increase at higher densities. Such a 
non- linear behavior has been observed, e.g., in ref. [22]. To show that this is 
a generic feature of such a type of models we calculate df(n)/dn, where f(n) 

5 In our case a difference appears for the u-meson and nucleon masses. 



$jv = $ CT = $ w = $ p = $(/) . 



(43) 



^ = $ 2 (/), 



i = 



(44) 
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is a solution of eq. (28) for isospin symmetrical nuclear matter with the total 
density n = n p + n n . From (28) we have 



df 



dn 



d 2 E N [n-f] 
dndf 



d 2 E N [n-f] 
dP 



f(n) 



(45) 



where the denominator d 2 E^jdf 2 should be always positive to assure the 
stability of the mean-field solution (28). The nominator, 



changes the sign at some value of n, provided drj^/df and d^w/df have 
the same sign. Indeed at small densities the second term dominates, whereas 
with the density increase the first term begins to win. Thus, if / increases at 
small densities (hadron masses drop), at larger densities / begins to decrease, 
resulting in an increase of hadronic masses. This effect can be enhanced in 
the isospin asymmetric nuclear matter due to the rho-meson contribution, if 
drjp/df has the same sign as drj^/df and d&N/df . The second feature of 
the energy-density functional (17) with scaling factors (44) is that it does not 
guarantee the continuity of the density dependence of the scalar field /. The 
equation of motion (28) for / may have several solutions, which can appear or 
disappear with a change of the density, leading to jumps of the f(n) function. 
Third, the scaling of the mass of the scalar cr-meson field included in a Dirac- 
Brueckner-Hartree-Fock calculation of nuclear matter leads to over-binding 
of the nuclear matter, as has been found in ref. [37]. Fourth, in ref. [38] it 
was noticed that increasing the vector repulsive potential, due to the factor 
1/Vu>, would produce too strong nucleon flow in heavy-ion collisions. Such 
a flow would be in odd with experimental data. This observation stimulated 
authors of ref. [20] to propose the scaling of the vector-meson coupling constant 
Xuj ~ &u>, which leads to r]^ ~ 1 . In the scalar sector the problems can be 
solved if we put r) a ~ 1 and choose the potential U so that there is only one 
solution of the equation of motion for /. 

Information about density dependence of the symmetry energy a sym (n), cf. 
(32), is rather controversial. The microscopic calculation within the relativis- 
ts Dirac-Brueckner-Hartree-Fock approach [35] shows that a sym (n) increases 
almost linearly with the density. This would require rj p ~ 1 in (32) . Ref. [25] 
discussing heavy ion collision properties uses a stiffer density dependence of 
a sym(^) (Vp decreases with a density increase). On an other hand, the HHJ EoS 
[16] fitting the microscopic model of Urbana-Argonne (Al8+5v +UIX*)[2] pro- 
duces a sym (n) oc n os that corresponds to increasing rj p ~ (n/n ) 0A . Note that 
the latter case allows to solve the problem of the low DU threshold density. 



d 2 E[n;f] 
dndf 



f(n) 



m 2 N rft df df y/m% ' 



C 2 n d Vw d<S> N m 2 N <S> N (f) 



(46) 
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3.2 RMF models with a universal scaling of hadron masses and coupling 
constants 



3.2.1 Construction of the model 

Following the conjecture of ref. [20] and argumentation in the previous sub- 
section we study now the case of a universal scaling of hadron masses and 
coupling constants, 

§ N = $ CT = $ w = $ p = x* = Xu = X P = $(/) , (47) 
which implies 



Va = Vu =V P = 1 ■ (48) 

Note that according to a hypothesis of a "vector manifestation" [39] the vector 
meson mass must go to zero in proportion to the corresponding coupling. At 
the saturation nuclear matter density refs. [20,40] give $(no) — 0.78 according 
to their analysis of the gyro-magnetic ratio that agrees with our input value 
(37), m* N (n )/m N ~ 0.8. 

Note also that for the case ^ ^ 1 we have a non-trivial relation between / 
field and an original a field (22). The W and ZM models (in both of them 
U — 0) with the universal scaling (47,48) we will denote as the W(u) and 
ZM(u) models, respectively, and the MW model with the universal scaling we 
will denote MW(u). For the W(u) and MW(u) models we then have 



*(/) = !-/ - f= g -f^r - *M = Y-W (49) 
For the ZM(u) model the correspondence is a more complex: 



$(/) = (l + /)- 1 — / 



/I 9aO_ _ 1 

4 m N 2 

-l 



^Hi\ +9 -^ + \) • (50) 



We note the obvious equivalence between the original Walecka model (W) 
and the W(u) model, and also between the original ZM and ZM(u) models 
for constant mean fields. We will denote models with a non-universal scaling 
as (nu) models. From (49) we see that W(u) model is in turn equivalent to 
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the ZM(nu) model with $ CT (cr) = &n(v) = 1/(1 + g a cr/mjv) but x„ — \ and 
^u;,p = — 1- Analogously we find equivalence of the ZM(u) model and the 
W(nu) model with ® atU1:P = Xuj, p = 1 and Xa(&) = ®n(<?) = 1 - 9aCr/m N . 

As it follows from (10,11,12,13,14), and (47) in the absence of the charged p 
meson mean fields and the potential U, with the help of the replacements 



we can equivalently transform a RMF model with a universal scaling of meson 
masses and coupling constants (r) a ,u>,p = 1) to a model without any scaling 
of meson masses and coupling constants. The absence of the non-Abelian 
interaction of p meson fields and the potential U is very important in order to 
arrive at such a conclusion, since the transition (51) would affect the non-linear 
sectors of the model. 

For the MW model with U given by (36) we can construct an equivalent 
MW(u) model with the universal scaling of meson masses and coupling con- 
stants (47) for $ = 1 — / and the same form of U written in terms of an / 
field. In terms of the a field the potential of the new MW(u) model looks more 
cumbersome: 



Note that potentials U(a) in the MW and MW(u) models, cf. (36,52), look 
quite differently in terms of the a field although in terms of the / fields they 
are the same. We used the same letter /, since the mean values of those / 
fields, determined from the very same equation of motion, are identical. This 
illustrates that, if thermodynamic characteristics look the same in terms of 
the variables found in the minimization procedure, then the RMF models are 
equivalent. 

As an illustration we also exploit a simplified ansatz proposed in ref. [20] 
(further the SBMR model) for the isospin symmetrical matter, where the 
Brown-Rho scaling was included into RMF EoS. Assume that the energy 
density of the model takes the form that we have used above with 



T)u = $ 2 (n;z)/$ 2 (n;2/), r) a = $ 2 (n; y) , $ N (n, f) = $(n; y) - f , 
$(n;y) = l/(l + yp|(n)/(260 MeV) 3 ) pj»(n) = (^+^)/2. (53) 

Here / = her/in^, h is the coupling constant and a = (ft in notations of 
[20], eq. (32). Parameters of the model are m N = 939 MeV, y = 0.28, z = 



H ^ € , where f = a, u, p , 



(51) 




b m N (g a a) 3 c (g a (r) 4 



(52) 
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Fig. 2. The effective nucleoli mass (left panel) and the energy per particle 
(right panel) of the isospin symmetrical nuclear matter within the MW(u) (in- 
put-parameter set (29,37,38)), ZM(u) (input-parameter set (29,40)) and SBMR 
model (53). 

0.31, Cl = 332.27 and C 2 a = 50.546, which correspond to n = 0.151 fin" 3 , 
e B = -16.1 MeV, K = 259.6 MeV, and m* N (n ) = 0.675 m N . Here is 
very closed to unit and r] a depends on n rather than on /. The pressure and 
chemical potentials of partilces are related to the energy following eqs. (15). 
Note that owing to an explicit dependence of the energy on the nucleon density 
expressions for the neutron and proton chemical potentials deviate from those 
given by eqs. (A. 10) and (A. 11). 

In order to consider asymmetrical matter within the same framework we add 
the rho-meson contribution (27) with i] p = 1 and C 2 p = 100.08. The latter 
quantity follows from the value of a sym (n ) = 32 MeV in (29). 

In Fig. 2 we show the effective nucleon mass and the EoS for the isospin 
symmetrical nuclear matter calculated for three models of EoS: MW(u) (input- 
parameter set (29,37,38)), ZM(u) (input-parameter set (29,40)) and SBMR, 
(53). We see that these three models cover a broad range of possible density 
dependences of the effective mass and the energy per particle. Compared to the 
RMF models the SBMR model produces very rapid decrease of the effective 
nucleon mass and a stiff EoS (at a reasonable value of the compressibility for 
n = no). 

3.2.2 Neutron star composition 

The composition of the NS matter (in absence of mixed phases and a charged 
condensate) is determined in the standard manner by the local charge-neutrality 
condition and the conditions of the equilibrium with respect to weak reactions 
e <-> p, and n <-> p + e . The latter ones impose chemical equilibrium conditions 
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n/n n/n Q 

Fig. 3. The same as in Fig. 2 but for the neutron star matter, 
for the corresponding chemical potentials 

« e = /V = « n - « p , Hi = Q n , i = e,fi,n,p. (54) 

Here E tot is the full energy-density. For the energy density (17) conditions (54) 
and the local charge-neutrality condition result in the system of equations 



He 



(" 2 



2\ - 

•e m e) 2 



3 7T 2 



id;) (>u< 



2m 2 N V P (f) 



3tt 2 



m 



(55) 



with / = f(n) being a solution of the equation of motion (28) for the scalar 
field. Eqs. (55) are solved with respect to n p for a given total density n = 
n p + n p . In Fig. 3 we show the effective nucleon mass and the EoS for the 
NS matter calculated for MW(u) (input-parameter set (29, 37,38)), ZM(u) 
(input-parameter set (29,40), and SBMR models. Compared to the isospin 
symmetrical case (Fig. 2) the EoS of the NS matter becomes substantially 
stiffer for all three models. 



In Fig. 4 (left and middle panels) we show the composition of the NS mat- 
ter and the electron chemical potential, as they follow from the solution of 
equations (54), (55) for the three specified above models (MW(u), ZM(u), 
and SBMR). We see that all three models produce the very similar NS com- 
positions and He ■ The proton concentration Y p = n p / n is rather high so that 
the DU process becomes operative at n — 2.13 no, 2.07 no, and 1.72 n for 
the ZM(u), MW(u), and SBMR models, respectively. The corresponding NS 
masses are M = 0.88, 0.99, and 0.97 M Q . Thus, yielding so low DU threshold 
densities these models cannot adequately describe the NS cooling, cf. [10]. This 
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Fig. 4. The proton and neutron concentrations (left panel), and the electron chemical 
potential (middle panel) in the neutron star matter calculated for the same models 
as in Figs 2, 3. The right panel shows the neutron star mass as a function of the 
central density. Horizontal cuts indicate the DU thresholds. 
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Fig. 5. The threshold density for the DU process as a function of the m* N (no) for 
the neutron star matter, for the different MW(u) models, the input set (29), at two 
values of K and three values a s 



is a serious drawback. In the right panel of Fig. 4 we show the NS masses (in 
units of solar mass M = 2 x 10 33 g) resulting from the integration of the 
differential Oppenheimer-Volkoff equations [1] for a given central density of 
the star. We stopped the integration at the inner crust-core boundary, be- 
cause the EoS following (17) cannot be used any longer. In different models 
the value of this boundary density varies in the interval n cmst — (0.5 -j- 0.8) Uq. 
In our calculations we settle n crust ~ 0.8 n . Since the crust contributes only 
a little to the NS mass (yielding an additional mass < 0.06 M ), we do not 
supplement our calculations by subsequent calculation of the crust. We see 
that the ZM(u) RMF model has an extra shortcoming that the limiting NS 
mass is rather low. Opposite, the SBMR model shows an example of a very 
stiff EoS (M lim ~ 2.8 M ). 

It is interesting to study the dependence of a DU threshold density and a 
limiting NS mass on the input parameters of a RMF model. We will use 
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Fig. 6. The counter plot of the limiting masses of neutron stars (counter labels in 
Mq), for the MW(u) models, the input set (29), for different values of the input 
parameters K and m* N (no). 



the MW(u) keeping the basic input parameters (29) and varying parameters 
m* N (n ) and K. Fig. 5 shows threshold densities for the DU process as a 
function of the magnitude of the effective nucleon mass at hq for two values 
of the compressibility K = 200 and 300 MeV and for three values of the 
symmetry energy coefficient a sym (n ) = 28, 32, 36 MeV. We see that 
increases with increase of K, and m* N (n ) and decreases with increase of a sym . 
In all cases is more sensitive to the changes of a sym and m* N (n ), than to 
the change of the compressibility. We see that all RMF models yield n^t < 
2.7 no for reasonable values of parameters. Hence all these models will meet 
difficulties in applications to a NS cooling problem, cf. NS masses in Fig. 4. As 
we have mentioned, models with M(n™ t ) < 1.3 M & meet serious problems, 
since the latter value is smaller than the averaged value of the NS mass, 
M ~ 1.35 ± 0.04 M Q , and models producing M(n™ t ) < 1 M should be 
completely excluded. 



The dependence of the limiting NS mass on the input parameters is illustrated 
in Fig. 6. The limiting mass is a sharp function of m* N (no) and a rather smooth 
function of K. We see that the limiting mass > (2.2 2.3)M Q can be reached 
either for models exploiting unrealistic values of the compressibility modulus 
K > 300 MeV, or for m* N (n ) < 0.7rriN. Therefore if one observed a NS of a 
mass M > 2.2 M , see [12], this would mean that many RMF models with- 
out any scaling or with a universal scaling of masses and coupling constants 
would have a problem. One of the possibility to avoid the problem is to use 
m* N (n ) < 0.7mjv- However, in the latter case one may meet a problem with 
description of the atomic nucleus data. Besides then either the compressibility 
is unrealistically high or c < 0, as it follows from Fig. 1. 
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3.3 RMF models with a non-universal scaling for masses and couplings 



3.3.1 Solving problems with DU reactions and limiting NS mass 

Obviously the universal scaling law may hold only approximately. In this sec- 
tion we demonstrate that allowing some departure from the universal scaling 
of the hadron masses and the coupling constants (47) we are able to obtain 
a somewhat stiffer EoS and shift the threshold density for the DU process 
to higher densities. We will consider the MW(nu) models with r\ c = 1 but 
Vu(f) < 1) being decreasing with n for n > n . This behavior enhances the 
repulsion in the potential at large densities. In absence of negatively charged 
bosons to increase the critical density for the DU processes we should reduce 
the proton concentration. This can be achieved by choosing rj p (f) > 1, being 
increasing with n for n > no, that results in a suppression of the rho- meson 
contribution to the EoS. 

As we have argued in section 3.1, if rj' w (f) < 0, then a decrease of a hadron 
mass will stop at some density and turn to an increase. However we can choose 
the scaling to be weak enough to shift a turning point to a high density not 
reached in the center of the heaviest NS (of the limiting mass). In the case 
of the NS matter we can exploit the fact that the first term in (46) respon- 
sible for the turn-over of the hadron mass scaling can be partially compen- 
sated by the analogous term from the rho-meson contribution, if r)'(f) > 0. 
For the pure neutron matter the exact cancellation condition takes the form 
^C^vLif) Vp(f) = ~^pVp(f)Vw(f)- An appropriate behavior of the scaling 
factors is obtained, e.g., with 



v* = i> = ' ^ (/) = 77V-7^T m w (56) 

where fo — 1 — m* N {no) /m^ • Since n^/o) = Vp(fo) = 1, the values of param- 
eters C w and C p do not change. Other parameters related to the scalar field 
change just slightly. 

We fit now the parameter z in order to increase the limiting NS mass Mn m 
and to push the threshold density to higher densities. The optimal value 
is found to be z = 0.65 . We use the same parameter set (29), as before. Had 
we taken the same values of the compressibility modulus K = 270 MeV and 
m* N (n )/m^ = 0.8, as we have used above, we would obtain negative value 
of c. The "c = 0" border shown in Fig. 1 shifts slightly up, if we apply the 
scaling (56). In order to get positive c we take slightly different values for the 
effective mass and the compressibility modulus compared to those in (38). We 
use 
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Fig. 7. The line shows the compressibility K of the isospin symmetrical nuclear mat- 
ter vs. the effective nucleon mass m* N {no) corresponding to c = for the MW(nu) 
model with the scaling (56) and z = 0.65 for the input set (29). 

m* N (n )/m N = 0.805, K = 275 MeV. (57) 
Then we obtain 



for MW(nu) z = 0.65 : 

Cl = 87.600 , Cl = 100.64 , C 2 = 179.56 , 

b = 7.7346 x 10" 3 , c = 3.4462 x 10~ 4 , (58) 

and 

for MW(u) z = : 

Cl = 87.600 , C 2 p = 100.64 , C 2 a = 184.36 , 

b= 5.5387 x 10~ 3 , c = 2.2976 x 10~ 2 . (59) 

In Fig. 7 we show a new c = border on the (K, m^)-plane for the MW(nu) 
models using the input set (29) and the scaling functions (56) for z = 0.65 . We 
draw the same conclusion as in the case of MW(u) models. In the framework of 
MW(nu) models (with z = 0.65) keeping c > we can't significantly decrease 
the effective nucleon mass m* N (no) not increasing the compressibility beyond 
realistic values. 

A very small value of the c-parameter in (58) that we have obtained, motivates 
a conjecture that the number of free parameters of the model can be still 
reduced. If we did not care about the global stability of the ground state we 
would put c = in (58). Taking care of this problem we could construct a 
model assuming 6 = 0, but c > 0. Exploring ^ ^ 1 we could put b = c = 0. 
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Fig. 8. The effective nucleon mass (left panel), and the energy per particle (right 
panel) of the isospin symmetrical nuclear matter (ISM) and the neutron star matter 
(NSM) within the MW(nu) model with the scaling (56), z = 0.65, input-parameter 
set (29,57,58). The dashed lines correspond to the MW(u) model [z = 0) with the 
input-parameter set (29,57,59). 

In this paper we will not consider these possibilities of further modernizations 
of our models. 

Fig. 8 (left) demonstrates the density dependence of the effective nucleon mass 
for the isospin symmetrical nuclear matter (ISM) and for the NS matter (NSM) 
for the given MW(nu) model (solid lines: z = 0.65, see (29,56,57,58)), and for 
the MW(u) model (dashed lines: z — 0, (29,56,57,59)). Dashed curves for 
ISM and NSM cases cannot be distinguished. We find that masses are slightly 
smaller for the MW(nu) model. We observe that for z = 0.65 the effective 
nucleon mass is monotonously decreasing in the symmetrical nuclear matter 
up to the density of 7 no- The effective nucleon mass becomes to increase only 
at higher densities. A dis-balance of the universal scaling is just minor. Right 
panel in Fig. 8 shows energies per particle. They are higher in case of the 
MW(nu) model (for n > 2 no). 

In Fig. 9 (left) we show how the given MW(nu) model (the input-parameter 
set (29,56,57,58)), may help to cure a problem of the low critical density for 
the DU process. For the MW(nu) model the DU process becomes forbidden 
up to 4.0 uq, whereas for the given MW(u) model (the input-parameter set 
(29,56,57,59)) we have ~ 2.08 no- The middle panel demonstrates the 
behavior of the electron chemical potential. The right panel shows the NS 
mass as a function of central density. With the MW(nu) model (z = 0.65) the 
limiting NS mass is increased by about 5%. We see that due to a stiffening 
of the EoS we may now satisfy the mentioned limit M > 2M . If we used 
another parameter choice, e.g., with a smaller effective nucleon mass at n , 
we could construct an even stiffer EoS. The critical value of the NS mass for 
the occurrence of the DU process is M(n™ ~ 4.0 no) — 1-69 M . This value 
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Fig. 9. The proton and neutron concentrations (left panel), and the electron chemical 
potential (middle panel) in the neutron star matter for the very same models as in 
Fig. 8. The right panel shows the mass of the neutron star as a function of the 
central density. The DU thresholds are indicated by cuts. 

is essentially higher than the averaged value of the NS mass measured in the 
NS binaries. Thus we removed the problem with the low critical density for 
the DU process. 



3.3.2 Fit to A18 + 5v + UIX* and EE J EoS 

Now let us show the efficiency of the MW(nu) models on another example. 
In ref. [16] a parameterization of the nuclear EoS was proposed (we called it 
the HHJ EoS), which is a good fit to the realistic EoS of the Urbana-Argonne 
group [2] (Al8+<5t> +UIX*) up to 4n . For larger densities the Urbana-Argonne 
EoS based on non-relativistic interaction produces a supersonic velocity. The 
HHJ EoS treats this problem, smoothly incorporating the causality constraint. 
In this parameterization the energy density is given by 



E[n p ,n n ] = (n p + n n ) 
u = (n p + n n )/n , e B 



e B u 



2.2 -u 
1 + 0.2w 



+ a. 



sym 



U 



0.6 \ U P 



45.8 MeV a sym = 32 MeV . (60) 



In ref. [11] it was demonstrated that the HHJ EoS can be successfully fitted 
by the MW model, which uses as an input set 



e B = -15.8 MeV , K = 250 MeV , m* N (n ) = 0.8 m N , 
a sym = 28 MeV . 

(61) 

The corresponding parameter set 
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Fig. 10. Left panel: the energy per particle of the isospin symmetrical matter 
(ISM) and the neutron star matter (NSM), as a function of the nucleon density. 
Right panel: the pressure of the neutron star matter. Solid lines show calculations 
done within the MW(nu) model with the scaling (63) and z = 2.9 with the in- 
put-parameter set (61,62). Dashed lines correspond to the MW(u) model (z=0) for 
the very same input-parameter set. Dash-dotted lines are for the HHJ EoS (60). 



Cl = 195.10, Cl 
b = 8.6964 x 1(T 3 . 



90.911, C 2 = 77.276 



c = 8.1411 x 10 



-3 



(62) 



is slightly different from the original one used in ref. [11], where the nucleon 
mass was taken to be 938.918 MeV instead of 938 MeV, that we adopted here. 
The MW(u) model with such parameters gives approximately the same depen- 
dence of the NS mass on the central density as the HHJ EoS [10]. However, as 
other RMF models, this MW(u) model yields a low threshold density for DU 
processes, n™ t ~ 2.6no (corresponding to the NS mass M ~ 1.1M Q ), whereas 
for the HHJ EoS ~ 5.2no (M ~ 1.8M Q ). This shortage can be corrected 
if we allow for a non-universal scaling of the rho meson mass and the coupling 
constant, which results in 



•n* = vu> = i , v P (f) = • ( 63 ) 

For z = 2.9 the scaling ratio $ P /Xp changes from unity at n to ~ 1.4 at 
7no- Thus a dis-balance of the universal scaling that we have introduced is 
not too strong. The threshold density for the DU processes is now shifted to 
n crit = 5.2 n , in accordance with that for the HHJ EoS. 

Fig. 10 (left panel) demonstrates the energy per particle of the isospin sym- 
metrical matter and of the neutron star matter (NSM) calculated within the 
MW(nu) model with the scaling (63) for z = 2.9 and for the input-parameter 
set (61,62) in comparison with the values calculated for the MW(u) model 
(z = 0), with the same input-parameter set, (61,62), and for the HHJ EoS. 
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Fig. 11. The proton and neutron concentrations (left panel), the electron chemical 
potential in the neutron star matter (middle panel) and the mass of the neutron 
star as a function of the central density for the very same models as in Fig. 10. 

Dashed and solid curves are not distinguishable. The right panel shows the 
pressure, as a function of the density for the NSM. We see that in both cases 
(ISM and NSM) the thermodynamic characteristics of the given MW(nu) and 
MW(u) models are very closed to the values calculated for the HHJ EoS. 

In Fig. 11 (left and middle) we illustrate how the scaling (63) allows to change 
the composition of the NS matter. We see that for the given MW(nu) model 
(input-parameter set (61,62)) the proton concentration and the electron chem- 
ical potential are strongly reduced compared to the corresponding MW(u) 
model (the same input-parameter set) and become equal to those for the HHJ 
parameterization. Right panel shows NS masses as functions of the central 
densities for the given three models. They are very close to each other. 

Concluding this section, we have shown an example of the MW(nu) model 
that allows to increase the limiting NS mass and to shift the DU threshold to 
higher densities. We have also constructed another MW(nu) model that can 
be used as a RMF model equivalent to the HHJ parameterization, thus also 
demonstrating an appropriate fit to the A18 + 8v + UIX* EoS. 



4 Second order phase transition to p condensate in NS 

Now we turn to the discussion of the isospin asymmetrical matter where, as 
we will show, mean fields of charged p mesons may have finite values. As we 
have mentioned, usually only one component, p^ is assumed to be non-zero. A 
new ansatz for the p meson fields was proposed in [26]. Besides the traditional 

(3) 

pi / field it includes also spatial components of charged p meson fields, 
pf = (pi ± ipf')/^2 7^ 0, i = 1,2,3. As it was argued in ref. [26], the 
presence of the fields p[ 3 \ p^ results in an increase of the system energy. 
Therefore, these fields can be put equal to zero. It was also found that the 
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condition (pfpj — pT Pj) — minimizes the energy. This implies that the 
ratio pf / pi is constant, independent of the spatial index i. Then we may take 
p~ = di p c and pi = di pi, where a — {a^} is the spatial unit vector, and p c is 
a complex amplitude of a charged p-meson field. 

(3) 

In terms of the p and p c fields the corresponding contribution to the density 
of the thermodynamic potential (effective energy) (21) becomes 



1 1 2 

^K. n p! /> Po\pc] = ^9 P X P (n P - n n ) pf - - (p{, 3) ) m 2 p <& 2 p 



(9 P X 



o {3) 



Pchj' 



m 2 p $1 



I Pc 



(64) 



Variation of the Lagrangian/ (thermodynamic potential) with respect to the 

(3) 

fields Pq and p i yields equations of motion 



(9 P Xppf -P p ch ) 2 -m 2 p ^ p c = 0, 

2 xJZ (3) . n ' ( I (3) p \ I |2 1 / \ 

m P $ P Po + 2 #p Xp (#pX p Po - P-ch) \Pc\ = ~9 P Xp K - n n ) 



(65) 



This system of equations has two solutions. The first solution is the traditional 
one, cf. (26), (27), 



Po 2^$^ (np ~ nn) ' Pc = °' 



,(3) _ 1 9p_ Xp_ 

P P 

yielding the contribution to the density of the thermodynamic potential 



(66) 



n(i)r n n-f} = C ^ Un - np)2 



(67) 



The second solution of (65) (found previously in [26], neglecting p p ch1 and 
without scaling of couplings) is as follows 



(3) Pch m p § p 2 
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b p 11 l N 'lp 

~cf 



X p 



v m P®P 



(68) 



It exists only, if \n n — n p \ > n p . The density of the thermodynamic potential 
corresponding to this solution takes the form 
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K, n pi f] = fi« [n n , n p ; /] - — ^— (|n n - n p | - n") 2 . (69) 

Hence, for n > n£ ,n , where the value n£ ,n follows from the solution of the 
equation 



n p = \n n - n p \ , (70) 

the second solution (69) becomes energetically favorable compared to the tra- 
ditional solution (66). The density n£' n which can be written as 



is the critical density with respect to a second-order phase transition. In prin- 
ciple a transition to a new solution could be of a first order with an abrupt 
change of the isospin composition and the density. It would then happen at 
a smaller critical density n = < n^ ,u , cf. [26]. We studied this possibil- 
ity (see below) and found out that in the presence of electrons the first-order 
phase transition is not realized. Thus we conclude that in the frameworks of 
the RMF models that we exploit here and with taking into account of electrons 
the charged p-meson condensation occurs by the second-order phase transi- 
tion. Note that the new solution (69) induces a spontaneous breaking of the 
spatial symmetry of the system, since we fixed only one value p c , whereas the 
choice of components Oj is not fixed. 

The charged p-meson chemical potential is found from the chemical equilib- 
rium conditions. Considering the reaction e + n <-> p~ + n we find 



/4 = /V (72) 

The number of protons, electrons, muons and p~ is governed by the local 
charge neutrality condition 



n p -n e -n, + n p ch = 0, n p ch = 2 (g pX ' p p^ - tf h ) \p c \ 2 < . (73) 
In the phase with p c ^ from (17), (69) we find 



1 

x (yfi& $jU/) + Pin ~ y/™ 2 N $2 v(/) + Pl P + ™ P ) > (74) 
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Fig. 12. Densities n p from (69) and n n — n p calculated for the neutron star matter 
within the MW(u) (input-parameter set (29,37,38)), ZM(u) (input-parameter set 
(29,40)) and SBMR models. Left panel shows the results for the choice x' = 1 
and right panel, for x' p = Crossing points of curves n p with the curve n n — n p 
correspond to n = n£' , cf. eq. (70). 
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The main uncertainty of our consideration is the unknown modification of the 
non-Abelian p-meson self-interaction in nuclear medium, which in our scheme 
is encoded in the scaling parameter x' p - As we have mentioned, we see no 
reason for the scaling of ^' thus taking x' p — 1- However for the completeness 
of the consideration and in order to investigate the sensitivity of our results 
to the parameter variation we consider also the case x' P — Xp — ^p- Moreover 
note that some QCD motivated studies, cf. [27], use smaller values for the 
p meson coupling constants g p = g ~ wip/F^ ~ 5.8, F w ~ 132 MeV, than 
those follow from the fits in the RMF models. In this respect the second choice 
X f p — xp — &p effectively demonstrates what could be if g p were decreased. 

The quantity n p (n) given by (69) and depicted in Fig. 12 for 

three types of models MW(u) (input-parameter set (29,37,38)), ZM(u) (input- 
parameter set (29,40)), and SBMR discussed above. We use here m p = 770 MeV 
The actual value of the critical density, n p ' 11 , for the second order phase tran- 
sition to the charged p-meson condensate state is given by the solution (71) 
of eq. (70). We see that the condensation is possible at n p ' u = 3.2 n for 
the MW(u) model (would be 4.1 no for the model with the scaling of yJ, at 
n p,n = 4 5 Ho f or £ ne 2M(u) model (5.6 no for the model with the scaling of 
X p ), and at n p ' u = 1.7 n for the SBMR model (2.2 n for the model with the 
scaling of x p ). 
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Fig. 13. Densities n p from (69) and n n — n p calculated for the neutron star matter 
within the MW(nu) model (63), z = 2.9 (the input-parameter set (61); (62)) 
are shown by solid lines. Dashed lines are for z = and the same choice of other 
parameters. Fat dots mark the critical densities n£' n . Curves 1 correspond to \ p = 1) 
whereas 2, to \ p = 

In Fig. 13 we show the same, as in Fig. 12, but for the MW(nu) model with 
z = 2.9 specified by (61,62,63). As the consequence of i] p > 1, the MW(nu) 
model produces a higher critical density for the charged p meson condensation, 
n p,n ^ 3_g Uo f or y _ i anc [ n p,n ~ 5.1 no for x P — & P - This is basically due to 
the smaller value of the asymmetry energy in (61) and therefore, the smaller 
value of the coupling constant C p . We point out that in a most realistic cases 
of the MW(u) and MW(nu) models the novel charged p meson condensation 
appears already for sufficiently low densities n ~ (3 -j- 4) n®. As the price for 
the increase of the threshold density for the DU process, the critical density for 
the condensation increased in case of this MW(nu) model. The MW(nu) model 
(56) for z = 0.65, using the input-parameter set (29,57,58) yields n£' n ~ 3.4 n 
for Xp — 1 an d ^c' H — 4-6 for x p — 

The critical density, n£' n depends on the particular values of the effective 
nucleon mass and compressibility at saturation, which are used to constrain 
parameters of the RMF model. We study these dependences on the example of 
the MW(u) model, which allows for independent variation of the input param- 
eters K and m* N (no)/mN- We calculate the critical density function 
of the input parameter m* N (no) / for K = 200 MeV and 300 MeV for two 
choices of the x! scaling mentioned above. Results are presented in Fig. 14. We 
see that n£' n is weakly dependent on K but varies strongly with the effective 
nucleon mass, especially for m* N (n )/m N > 0.75. For 0.6 < m* N (n )/m N < 0.8 
the critical density varies within the interval Uq < n£ ,n < 3n for x p — 1- 

Fig. 15 demonstrates the energy per particle and the pressure of the NS matter 
as a function of the nucleon density without (dash curves) and with p~ con- 
densate (solid curves) for the case of the MW(u) model (the input-parameter 
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Fig. 14. The critical density of the charged p meson condensation, n£ ,n as a func- 
tion of the m* N {no) used as an input of the MW(u) model for two values of the 
compressibility modulus K and three values of a sym . Other input values are taken 
according to (29). The solid lines are calculated for the choice x' P = 1 an d the dashed 
ones for \' p = 




n/n n/n 

Fig. 15. The energy per particle (left panel) and the pressure (right panel) of the 
neutron star matter as a function of the nucleon density without (dash curves) and 
with (solid curves) p~ condensate for the MW(u) model (the input-parameter set 
(29), (37); (38)). The lines labeled by "1" are calculated for the choice x' P = 1 an d 
the lines labeled by "2" are for x' p = ^ P - 

set (29,37,38)) cf. Figs. 2, 3. Curve 1 is for x p — 1 an d curve 2, for \ p — & P - 
We observe no van der Waals behavior of the pressure typical for the first 
order phase transitions. The EoS is softened in the presence of the second 
order phase transition to the charged p~ condensate but quantitatively the 
softening effect is minor. 

In Fig. 16 we show the proton fraction Y p = n p /n, the electron chemical 
potential, and the NS mass for the NS matter as a function of the nucleon 
density without (dash curves) and with (solid curves) p~ condensate for the 
case of the MW(u) model (the input-parameter set (29,(37,38)), cf. Fig. 4. 
Curves 1 are for x p = 1 and curves 2 are for x p — $ P - We see that a part 
of electrons is now replaced to the condensate p~ mesons. Due to that the 
electron chemical potential decreases in presence of the condensate compared 
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Fig. 16. The proton concentration (left panel) and the electron chemical potential 
(middle panel) of the neutron star matter as a function of the nucleon density 
without (dash curves) and with (solid curves) p~ condensate for the MW(u) model 
(the input-parameter set (29,37,38)). The right panel shows the mass of the neutron 
star as a function of the central density. The lines labeled by " 1" are calculated for 
the choice x'p = 1 an d the lines labeled by "2" are for x'p = 



300 



500 




Fig. 17. The same as in Fig. 15 but for the MW(nu) model with (63), z = 2.9, the 
input-parameter set (61,62). 

to the case without the condensate. The proton fraction becomes smaller in the 
presence of the p~ condensate than would be in the absence of the condensate. 
However, as we obtain in the framework of the MW(u) model, DU processes 
once started at 2.07 no stay operative in the condensate phase up to rather high 
densities (up to 5.8 no for case 1 and up to the maximum central density for 
case 2). The neutron star masses change very little with the appearance of the 
condensate. 

Figure 17 shows the same as Figure 15 but for the MW(nu) model with (63), 
z = 2.9, the input-parameter set (61,62). Again we find no van der Waals 
behavior of the pressure and only weak softening of the EoS. The proton 
concentration, electron chemical potential and the NS masses for this model 
are shown in Fig. 18. The only qualitative difference is seen in the behavior 
of the proton fraction for the case x P = ^p- It increases with increase of the 
density in this particular case. However the triangle inequality for the DU 
processes is not fulfilled since a part of the negative charge is taken by the p~ . 
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Fig. 18. The same as in Fig. 16 but for the MW(nu) model with (63), z = 2.9, the 
input-parameter set (61,62). 




Fig. 20. The same as in Fig. 16 but for the MW(nu) model (56), z = 0.65, the 
input-parameter set (29,57,58). 

Thus in the case of the MW(nu) model, z = 2.9, DU processes do not occur 
in the phase with the condensate at all. 

Fig. 19 shows the same as Fig. 17 but for the MW(nu) model (56), z = 0.65, 
the input-parameter set (29,57,58). Fig. 20 shows the same as Fig. 18 for the 
given model. Qualitative behavior of the quantities for the given MW(nu) 
model (56) with z = 0.65 is the same as for the MW(nu) model (63) with 
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z = 2.9. Also in the case of the MW(nu) model, z = 0.65, DU processes are 
non-operative in the condensate phase for both cases x — 1 and x p = $ p . 

Concluding this section we found that a new phase of the dense nuclear matter 
with a charged p meson mean field might exist in the NS interiors already for 
not too high densities. For the most realistic cases we found n P ' n < 4 n . 



5 Neutrino emissivity in the DU processes on p condensate 



As we have demonstrated, the standard DU processes n — > peu, pe — > nv are 
forbidden in cases of the MW(nu) models (56,58) and (63,62) up to sufficiently 
high density (n > 4 n and n > 5.2 n for the case x p — !)• Therefore we 
consider these models as most realistic ones. However a new channel of the DU- 
like processes on the condensate charged p~ meson, e.g., np~ — > neu becomes 
operative for n > n^ 11 . Since the charged p~ condensate arises in the s-wave, 
the calculation of its contribution to the emissivity can be done similar to 
that for the kaon condensation processes, except additional requirement of the 
chiral symmetry used in ref. [46]. Note that ref. [46] did not take into account 
nucleon-nucleon correlations in vertices that yields an additional suppression 
factor < 10 _1 for all condensate processes, cf. [47,5,48]. Recovering the missing 
factor we easily estimate the neutrino emissivity 



e 



$ ~ 10 26 ^X? (^V (^) T 9 6 ^^r\ (75) 
p p \m N J \m N J m p <$> p m% cur • sec 



= 140 MeV, T 9 = T/10 9 K, T 4 = T 2 S T 2 W _ S is the nucleon-nucleon correlation 
factor, 



r s = r(/',e ~ n e ,k ~p Fn ), T(x,e,k) = [1 - 4C xA np (e,k)] \ 

T 2 w _ s = T 2 (f',e = q^p Fe ), (76) 

Co c± 0.77m~ 2 relates to the density of states at the Fermi surface for n = n , 
/' ~ 0.5 4- 1 is the Landau-Migdal parameter in isospin channel, cf. [49,5,50], 
A np is the proton - neutron hole loop (without spin degeneracy factor), see 
[48]. As we have mentioned, the presence of the correlation factor is crucially 
important. It may suppress the rate by a factor ~ 10 _1 4- 10 -2 . Another 
suppression factor arises from a small concentration n p ch in the new phase, 
n p ch < rip. 

In presence of the nucleon superfluidity the emissivity of the DU-like processes 
is additionally suppressed by a factor min[^ n „, ^ ~ e -Au/T f or j> < 
where A,, is the pairing gap in the corresponding nn or pp channel. According 
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to [51,10] the 3 P2 nn-pairing gap should be very small, less than 10 keV. If so, 
min[£ nn , £ pp ] ~ £ pp . Final rough estimation of the emissivity is 

e* ~ (10 24 - 10 25 ) Tl min[£ nn , U — ^— . (77) 

cur • sec 

This rate is typically either smaller or of the order of that expected for the 
charged pion condensate. However it is much smaller than that for the stan- 
dard DU processes, if the latter were permitted. On the other hand, for typical 
densities and temperatures T ~ 0.3-=-3- 10 8 K the emissivity on the charged p~ 
condensate is somewhat larger than the emissivity of the competing medium- 
modified Urea processes, as they have been computed in [52,10]. 

As has been shown in ref. [10] within HHJ EoS, the critical density of the 
charged pion condensate should not be less than (2.5 -r- 2.7)no- Otherwise 
even rather low mass NS would cool too fast in disagreement with the data. A 
similar statement can be done in case of the charged p~ condensation. Thus we 
see that the values n£ ,n ~ (3-j-4) no that we have obtained within our MW(u) 
and MW(nu) models do not contradict to the analysis of the NS cooling [10]. 



6 Conclusion 

In this paper we studied generalized relativistic mean field models that include 
in-medium changes of meson masses and coupling constants. We demonstrated 
equivalence of a certain class of the models. Practically we showed that the 
energy density for homogeneous matter depends only on three independent 
combinations of mass and coupling scaling functions. Often this equivalence 
is not realized and one considers some models as different, whereas in reality 
they are identical for the description of homogeneous matter on the mean field 
level. 

We showed that a naive modification of solely boson masses would lead to se- 
vere drawbacks, like discontinuous changes of mean fields and hadron masses 
with the density. To get rid off these problems the mass scaling should be es- 
sentially compensated by the scaling of the corresponding coupling constants. 

For linear RMF models, i.e., without a non-Abelian p interaction and a non- 
linear a potential U, we demonstrated equivalence of the models with universal 
scaling g* rrii/ giin* ~ 1 , % — a, u, p, and the models without any scaling of 
those masses and couplings. For non-linear RMF models (U ^ 0) the above 
equivalence would require a certain relation between corresponding non-linear 
potentials. 
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We have demonstrated the efficiency of the generalized relativistic mean field 
models based on the modified Walecka (MW) model with a non-universal 
scaling of meson masses and couplings (MW(nu)). It allows to enlarge the 
threshold density for the direct Urea process and produces a stiffer equation 
of state without any changes of the EoS near the saturation density. These 
possibilities are favored by "the nuclear medium cooling scenario" for the 
neutron star cooling, cf. [10] and by recent measurements of heavy neutron 
stars, cf. [12]. 

Our generalized mean field models with a non-universal scaling of meson 
masses and couplings might be useful in discussion of possibilities of different 
phase transitions in neutron star interiors, since they allow to compensate a 
softening of the equation of state due to a phase transition by a stiffening 
due to a stronger decrease of masses compared to a decrease of couplings. 
Especially the model might be helpful in discussion of hybrid stars. With the 
standard hadron models the quark matter either begins to appear already at 
rather low neutron star mass, or does not appear at all, cf . [45] . In the former 
case the limiting mass of the star is found to be rather low (typically < 1.8 M Q ). 
If one used a model assuming a smaller nucleon mass at the saturation, one 
could get a still higher value of the limiting neutron star mass. 

We also demonstrated an example of the relativistic mean field model that 
fits well the modern Urbana-Argonne equation of state (more precisely, the fit 
[16] that cures the causality problem), yielding the same threshold density for 
the direct Urea process. 

The MW(nu) models can be useful for the description of properties of finite 
nuclei. Here equivalence is absent due to the spatial inhomogeneity of the field 
profiles. In standard MW models one is forced to decrease the a meson mass 
in order to appropriately fit the proton density profile and the binding energy 
as a function of the atomic number, cf. [44]. In the framework of the MW(nu) 
models there is still a room to fit the scaling laws for effective masses and 
couplings. 

It would be interesting to apply the model to the description of heavy-ion 
collisions. Early naive studies of the pion yield in heavy-ion collisions required 
extremely hard EoS, with the compressibility K ~ 1000 MeV. Further it was 
demonstrated, cf. [5,53,54] and references therein, that the problem might 
be naturally solved within the standard modified Walecka model (MW) when 
one includes contributions from soft pion modes. These modes yield very small 
contribution to thermodynamical characteristics of the cold nuclear matter but 
produce a large contribution for the hot nuclear matter. The physical reason 
is obvious: soft modes are easily excited with the temperature increase. Due 
to a decrease of a, cu, p and nucleon masses related modes become to be softer. 
Thus our mean field models may give an appropriate basis to go beyond the 
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mean field level considering thermal excitations. 

Recent studies of the collective isospin flow [23] extracted the density depen- 
dence of the symmetry energy but did not introduce any astrophysical con- 
straints and a possible soft mode contribution. It is therefore a challenge to 
fit the corresponding data using the constructed above models with inclusion 
of thermal contribution from soft modes. 

If one introduces a rho- meson field as a non-Abelian gauge boson (cf. [27]) in 
a relativistic mean-field model which supports the decreasing p-meson mass, 
one should take care of the possibility of the charged rho meson condensation 
in neutron star matter. Following [26], we found that at density n > n£' n the 
charged p-meson mean field appears. We demonstrated that in the framework 
of our relativistic mean field models the novel phase arises by a second order 
phase transition. The critical density n£ ,n is rather low in the neutron star 
matter (typically ~ 3 4- 4 no), if one allows for a decrease of the rho meson 
mass. We also included the corresponding decrease of the effective nucleon- 
rho coupling, as well as the scaling of other meson masses and couplings. If 
one used a model assuming a smaller nucleon mass at the saturation, one 
could get a still smaller value of the critical density for the charged p conden- 
sation. Charged p condensation slightly softens the equation of state. In the 
neutron-enriched matter the charge of the p-meson condensate is negative. We 
evaluated the emissivity of the direct Urca-like processes on the condensate 
p~ meson. The resulting rate behaves similar to that for the n~ condensate. 
The analysis of the cooling data indicates that the value n£' n should be not 
too small, n p c ' 11 > 2.5 4 2.7 n . 

Our non-Abelian p meson field is in some sense analogous to the gluon field. 
Dropping of the effective p meson mass enforces a similarity. Effective nucleon 
and other meson masses also decrease with increase of the density. Decrease of 
the corresponding effective couplings motivates a future asymptotic freedom 
in the quark-gluon phase. All these features might be in favor of a smooth 
transition from the hadron degrees of freedom to the quark-gluon degrees 
of freedom. Up to now there is no appropriate description of a cross-over 
temperature region seen from the lattice results [55]. Formulated above models 
might be helpful in construction of such a description. 
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A Thermodynamic potential of RMF model in presence of the 
electric field 



So far we have disregarded the electromagnetic interactions. But only including 
the coordinate dependent electric field one may properly discuss structures, 
that may occur at first order phase transitions in charged systems. Indeed, a 
charged structure may exist, only if there is an inhomogeneous (not constant) 
electric potential. Although we did not find a first order phase transition in the 
framework of our RMF model EoS, this result might be still model dependent. 
Here we demonstrate how one can include electromagnetic interactions into 
the scheme. We consider the NS matter consisting of the neutrons, protons, 
electrons and muons. We incorporate mean meson fields including charged p 
meson field, p~ = a^p c , p~ = a ti p\ with a M = (0, a) and the electric field. 

The interaction with the static electric field eN 1 = (— V, 0) can be introduced 
according to the minimal coupling prescription via the gauge replacement 



d^ N -> (d^ + l -eA^ (1 + r 3 )) ^ N , d„% -> (fy - ieA^ , 

Ppv ~^P»v-e [n 3 x p J + e A v [n 3 x , (A.l) 

where {n 3 ) a = 5 a3 is the unit vector in the isospin space, a = 1, 2, 3 . The 
frequencies of the charged fields acquire the shift e — > p — V for negatively 
charged fields (e, p~ and p~ in our case) and e — > p + V for positively charged 
fields (for proton). 

The density of the total thermodynamic potential is as follows: 
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I -prV™? +P 2 -ni(»i-V) l = e,», (A.7) 
o 



for |/i e - y| > m M . 

With the help of the variation 



and 



(^j ttdr^j =0 i = n, p, e,/i (A. 8) 



(| fidr) =0, ^ = <r, u, pf, 3) , p c , K (A.9) 



we recover necessary equations of motion. In the body of the paper we used 
the same equations but for fields independent of coordinates. 

Chemical potentials of neutrons, protons, electrons and muons are found from 
equations of motion (A. 8): 

= $ 2 N (f) +p| ;n + g u x<J*o ~ ^gpXpPo* - ( A - 10 ) 



th + V = ^m 2 N $ 2 N (f)+pl p + guXuWo + \g P X P Po\ (A. 11) 



Hi 



~V = Jm[+pTy (A.12) 



Chemical potentials of proton, neutron, electron, muon and charged p are 
related to each other following equilibrium conditions in reactions n <-> p + e, 
n <-> p + fj,~,n <-*p + p~. Thus p n = p p + p e , p e = P P = pf ch . 

The local charge neutrality condition that we used in the paper body is now 
replaced by the Poisson equation for the electric potential 

AV = -47i(n p -n e - nil + n p ch ), (A.13) 
where 



[(V-Hw) 2 -mlJ*'* nmr 
n e , P = 9(\V - He,?] ~ m etll ) (A.14) 
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for electrons and muons and 



< h = (g P x' P p { o } +V- a&) \ Pc \ 2 < (A.15) 

for the charged p~ meson condensate. The local charge neutrality condition 
appears as the solution of the Poisson equation for the charged structure of the 
size much larger than the maximum value among Debye screening lengths of 
different species. For the case of Maxwell construction, solving (A. 13) one may 
describe the boundary layer between phases. With these findings one may also 
solve the problem of possible occurrence of the mixed phase structures. This 
discussion is in full analogy to the framework of ref. [44] describing mixed phase 
at the first order phase transition to the s-wave negative kaon condensate. 
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